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Abstract: If future colliders discover supersymmetric particles and probe their 
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We investigate how well the relic density can be predicted in minimal supergravity 
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tify the collider accuracy needed. Theoretical errors in the prediction are investigated 
in some detail. 



Keywords: park Matter, Beyond Standard Model, Cosmology. 



*address after 1st Oct: DAMTP, CMS, Wilberforce Road, Cambridge, CB3 OWA, UK. 



Contents 



|T|. Introduction 



I. mSUGRA post WMAP 
||. The set-up 



37l| The approach 

Deriving the accuracies through an iterative procedure 
Analysis in the mSUGRA approach 



3^ Analysis in the PmSUGRA approach 



Co- annihilation 

Characteristics of the co-annihilation sliver in mSUGRA 
Theoretical uncertainties in mSUGRA 



4.1 



4.2 



47^ Accuracies of the GUT scale mSUGRA parameters 



CI Accuracies on the relevant physical parameters 



Higgs Funnel 



5A\ Characteristics of the funnel region 

|5.2| Theoretical uncertainties in mSUGRA 

|5.3| Required accuracies on the mSUGRA parameters 

[5^ Pmsugra 

p. Focus point 



6.2 



0| The focus point landscape 
Theoretical uncertainties 
Accuracies in mSUGRA 
Accuracies in PmSUGRA 



6.3 



6.4 



1^. Summary and Conclusions 
1^. Inputs 



H 
u 

m 
m 

n 

21 
2| 
24 
26 

28 

28 
3G 
32 
33 

35 

|37| 



- 1 - 



1. Introduction 



One of the attractive features of the minimal supersymmetric standard model (MSSM) 
is that it provides a natural candidate for cold dark matter, the neutralino x? [Q]- 
With cosmology entering the era of precision measurements and the next colliders 
aiming at discovering and constraining supersymmetry some crucial cross breeding 
is emerging. Already, assuming the standard cosmology |@], the measurement of the 
relic density of dark matter has been used to put strong constraints on the supersym- 
metric model 0, ^ ^, 1^. For example, WMAP [^, which at lo constrains the relic 
density in the range .094 < fi/i^ < .129, effectively reduces the dimensionality of the 
MSSM parameter space by one. is the mass density in units of the critical density 
and h is the scaled Hubble constant. The accuracy of this constraint is expected to 
increase with future data from the PLANCK satellite, which should obtain precision 
on Qh^ at the 2% level g |TOl . 

Previously, much of the literature has used ^h? as well as other measurements 
in order to place constraints upon the MSSM [^, ^ ||, |Tl|, |T2|, [1^, [1^, sometimes 
giving predictions for the signatures expected at colliders and direct dark matter 
searches. In this paper, we look at the problem from the inverse perspective and 
examine what is required from collider data in order to get a precise prediction for 
n/i^, which could then be used to test the cosmology. For example, we can ask: is all 
of the dark matter due to the lightest neutralino, or is there some other component? 
Assuming that the neutralino constitutes all of the dark matter, the predicted fi/i^ 
value for neutralinos may test the cosmological assumptions that go into the deriva- 
tion of the relic density in the standard picture. A value predicted from colliders 
with enough accuracy that does not match with the WMAP number would indicate 
some non-standard cosmology. This would be highly exciting. 

In the standard picture, at early times, the would be dark matter relics are in thermal 
equilibrium with a number density n^q oc {MTY/'^exp{—M /T) for a non-relativistic 
particle of mass M at temperature T. Once the interaction rate of the particles falls 
below the rate of expansion of the Universe set by the Hubble parameter, the particles 
freeze-out with a number density determined by the so-called freeze-out temperature 
Ty, with Ty ~ Mj/25. Considering the exponential decrease of n^q the relic density 
therefore depends critically on the precise moment the annihilation rate equals the 
expansion rate of the universe. This interaction rate and its temperature/time evolu- 
tion is completely controlled by particle physics. On the other hand, in the standard 
approach one implicitly assumes that the post-inflation era is radiation dominated. 
There is no clear-cut evidence that between the freeze-out (at T ~ 10 — lOOGeV, 
say) and big bang nucleosynthesis BBN (at T ~ IMeV) the universe was radiation 
dominated. Radiation domination could start at a "reheat" temperature that may 
be quite low. If the expansion of the universe is modelled differently than in the stan- 
dard picture so that some other contributions or dynamics occur, the link between 
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particle physics and cosmology is no longer unique. As examples, one can mention 



a few explicit scenarios such as a low reheating temperature [|15|, |16|, scalar field 



kination|[T^, |T8], |T9|, scenarios with extra dimensions I^Tj, anisotropic cosmology 



3i, or scalar-tensor cosmologies E^. Entropy creation as well as non-thermal 



production of neutralinos can also change the picture. 

For the particle physics part, in order to predict the relic density Qh"^ in the 
MSSM one generally needs to know many of the numerous underlying parameters of 
the model, since that determines the variety of annihilation channels for neutralinos 
in the early universe. Indeed, in the standard approach Qh"^ oc 1/ < av >, where 
< av > is the thermally averaged cross section times the relative velocity of the LSP. 
How many parameters could be known and with which precision depend both on 
the theoretical model and the future terrestrial data, especially from the colliders. 
At future colliders such as the Large Hadron Collider (LHC) and a future linear 
collider facility (LC), supersymmetric sparticles are expected to be produced and 
measured if low energy supersjTumetry (SUSY) is present in nature p6|, It is 



hoped (though not guaranteed) that enough information |2^, ^] will be present in 
order to discriminate between various models of supersymmetry breaking. When we 
have identified a successful model of SUSY breaking, an accurate prediction of Qh"^ 
will allow us to test the cosmological assumptions that go into its prediction and to 
determine whether all of the dark matter is due to the neutralino. This identification 
of the underlying SUSY breaking scenario, although highly desirable, may not always 
be possible. However to predict the relic density it may also happen that one does not 
need to know all of the parameters but only those relevant parameters that are needed 
for the relic density calculation. An example of what we mean by an underlying 
SUSY breaking scenario is the minimal supergravity, mSUGRA^ model. Even when 
an underlying model has been identified, at the LHC for example, the properties of 
those particles which play a dominant role in the relic density calculation are not 
measured directly or precisely. The relevant parameters are then inferred through a 
complex theoretical framework that brings with it an associated uncertainty. 



One aim of this paper is to critically examine the theoretical uncertainties ^ 
in the calculation of the MSSM spectrum and their impact on a precise determina- 
tion of the relic density as well as deriving the accuracy needed on the fundamental 
parameters of the underlying model. We also take the model independent approach 
where we aim at only identifying the set of relevant parameters which are the most 
important parameters to measure in different regions of parameter space, and what 
empirical accuracy upon them will be required. We will refer throughout to the 
WMAP benchmark on accuracy: those that produce a 10% change in Qh"^ ('WMAP 
accuracy"). In part we also discuss what is needed to achieve the 2% level ("PLANCK 
accuracy"). Here, we will identify requirements on observables more generally (along 



Often referred to elsewhere as the CMSSM. 
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lines in WMAP-allowed parameter space) and although we leave speculations about 
the feasibility of meeting them to the conclusions, we discuss which additional col- 
lider observables are worth investigating to improve on the precision of the relevant 
parameters. In this article we carry this investigation in the mSUGRA framework 
with the supersymmetric spectrum provided by SOFTSUSYl . 8 . 7 The relic den- 
sity is computed with micrOMEGAsl . 3 P3[ which is interfaced with SOFTSUSYl . 8 . 7 



through the SUSY Les Houches Accord 



The emphasis upon corroborating accurate cosmological measurements with ter- 
restrial collider data, rather than cosmology constraining SUSY, has been addressed 
very recently, see for example ||35||. Ref. ||36| investigated how the variation of the in- 
put parameters affected the prediction of the relic density (and direct detection) but 
this was a pre-WMAP study in a scenario no longer viable and, more importantly, 
does not address the sensitive issues of poles and co-annihilations that we will study. 
Ref. examines how, at the LHC, one could determine some parameters entering 
the relic density calculation, while Ref. determines the required accuracy on the 
fundamental parameters of mSUGRA before translating this into an error on the relic 
density. Again both these studies were pre-WMAP studies and hence in a region of 
parameter space which is quite favourable for the LHC. Some recent experimental 
simulations covering part of the mSUGRA parameter space compatible with WMAP 
have been performed in ^ for the LC and in [|14| for the LHC. 

The paper is organised as follows. We will start in Section |^ by briefly describing 
the three regions in the mSUGRA parameter space that are still compatible with the 
WMAP constraint and which we will study in detail in the next sections. These 
are the f co-annihilation region, the Higgs funnel region and the focus point region. 
We felt that a dedicated Section ||, devoted to how we quantify the sensitivities and 
how we approach the problem, both purely within mSUGRA and in a more model 
independent way, was warranted. We investigate the co-annihilation region first in 
Section |. Results for the Higgs funnel region are shown in Section ^ In section ^ 
the focus point regime is studied. Conclusions follow in section |^. Input parameters 
for the analysis are detailed in Appendix |^. 



2. mSUGRA post WMAP 

mSUGRA is defined through only four parameters and a sign: Mi/2,mo, Aq i.e. the 
common gaugino mass, scalar mass and tri-linear coupling respectively (all defined 
at the GUT scale) as well as tan/5 the ratio of the vacuum expectation values in 
each of the Higgs doublet at the weak scale. The sign refers to the sign of the 
derived /i parameter, the Higgsino mass term. Just prior to WMAP, mSUGRA was 
compatible with the limit Qh"^ ~ 0.1 — 0.3 and other direct and indirect low-energy 
and collider data in a rather large region of parameter space called the "bulk region" . 
The SPSla mSUGRA benchmark point lies in this region. It was studied in 
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an ATLAS simulation with the assumption that nature follows mSUGRA[3S]. It 
was concluded that a precision of 3% on the relic density was possible in this case. 
After the WMAP constraint upon Qh"^ however, and given empirical lower bounds 
on sparticle masses, this large bulk region in rriQ — M1/2 is no longer viable, see 
for example |T^. SPSla is now ruled out by the WMAP relic density constraint. 
The point is that in mSUGRA, the neutralino LSP happens to be, for practically 
all cases, an almost pure bino that annihilates most efficiently into leptons through 
the right-handed sleptons because of their larger hypercharge assignment. However 
with the newest WMAP data this mechanism is not efficient enough. There remains 
then three favoured scenarios which all require some very specific accidental relations 
between some parameters at the electroweak scale. 



• f co-annihilation: 

In mSUGRA at small mo, there exists a region with almost degenerate f — x^. 
In this case the population of these two sparticles is almost the same, making 
the NLSP f thermally accessible. AM = — m^o is the mass difference which 
controls the ratio of the population of the two species through the Boltzmann 
factor exp(— AM/Tf) p^. It is then a very sensitive parameter that enters the 
calculation of the relic density. When co-annihilation takes place, through the 
participation of the fi in processes such as fiXi T7 or even fifi rf, the 
relic density can be brought down compared to the case of the bulk scenario. 



Higgs funnel region: 

A sudden increase in the usual annihilation mechanism to bring down the relic 
density can also occur if m^o is near a pole. Collider constraints on the LSP in 
mSUGRA allow the heavy "Higgs funnel" Hf: where XiXi ^ bb/rf, 
which occurs at large tan /3. 



• Focus region: 

Most of the time in mSUGRA, /i is rather large. However it may exceptionally 
happen that fi ~ Mi or even fi < Mi, in which case the annihilation is much 
more efficient through reactions such as XiXi ~^ WW/ZZ/Zh/ti. This occurs 
in the so-called focus point region ^ where tuq is very large. 



Our study will cover these three rather constrained regions in the mSUGRA 
parameter space in detail. For all the MSSM models considered here, the LSP is a 
neutralino. It is therefore useful to list its parameters through the neutralino mass 
matrix which will help understand some of our results. 

The neutralino mass matrix, in the bino, wino, "up" and "down" Higgsino basis 
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respectively is defined as 



( Ml -MzsmOw cos 13 sin 6*1^ sin /? \ 

M2 Mz cos 6w cos 13 —Mz cos 6w sin (3 

-M^ sin cos /3 cos 6'vi/ cos /3 — // 

y Mz sin 6*^/ sin [3 —Mz cos 6'vi/ sin [3 —fj, y 



(2.1) 



where Mi and M2 correspond to the weak scale bino and wino masses. Mz is the 
Z mass. This matrix is diagonalised with a matrix that we label N. In mSUGRA, 
with 6w the weak mixing angle one has 

5 

Ml = - tan^ 9wM2 ~ M2/2. (2.2) 

3 

This matrix receives radiative corrections but keeps the same form as the above 
tree-level structure. In our analysis these radiative corrections consist of an approxi- 
mation to the full one- loop result, based on Ref. All sparticle mixing is ignored 



in the correction terms, corrections proportional to the U{1) coupling are ignored, 
quark masses are set to zero and the squarks are approximated to be degenerate. 
Also, tree-level formulae for the Higgs are used in the loop corrections to the mass 
matrix. The tau Yukawa coupling is neglected. These approximations induce errors 
of order (a/47r)M|//i2 and (a/47r)M|/M|. Ma is the mass of the pseudo-scalar 
Higgs. 



3. The set-up 
3.1 The approach 

To "generate" the models compatible with the WMAP data we will fix Aq, tan (3 and 
the sign of but scan over the GUT value of the gaugino mass M1/2. The strong 
constraint from WMAP then means that mo can not be varied at will but, is to a 
large extent, determined by Mi/2- We therefore first seek a parameterisation of the 
mo-Mi/2 dependence which describes the three regions of the mo-Mi/2 plane which we 
are interested in. Each region will then be defined by a slope. Of course, in order to 
arrive at the full spectrum at the electroweak scale that is needed for the calculation 
of the various rates entering the relic density, one also needs to specify the standard 
model parameters, such as mb{mi,),mt etc. To derive the weak scale parameters we 
rely on SOFTSUSY. SOFTSUSY includes the necessary RGEs (Renormalisation Group 
Equations) and calculation of pole masses. It is well known that this procedure neces- 
sarily introduces some theoretical errors on the prediction of the weak scale physical 
parameters depending on the level of sophistication (loop order) and approximations 
that are implemented in the code. In a preliminary analysis we have compared dif- 
ferent codes for the spectrum evaluation and their impact on the prediction of the 
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relic density [50 1 . Here, for all three cases, we estimate the theoretical uncertainty 
solely with SOFTSUSY. For example, SOFTSUSY uses a renormalisation scale Msusy 
at which the electroweak symmetry breaking conditions are imposed, and at which 
most of the sparticle masses are calculated. By default, Msusy = is the 

renormalisation scale used at which it is hoped sensitivity to higher order corrections 
is small. In principle, changing this scale only has an effect on physical observables 
at a perturbative order higher than that used in SOFTSUSY. The present theoretical 
uncertainty coming from higher order effects is then estimated by re-calculating VlK^ 
for different values of 

^ = M'susY /Msusy- (3.1) 

Here, we take x = 0.5 — 2.0 to represent the scale variation. We will, depending on 
the situation, also study other uncertainties stemming from the different values of the 
SM parameters. Some of the theoretical uncertainty in SOFTSUSY could be reduced 
by using the full two-loop effective electroweak potential. Such a potential has been 
presented in Ref. ||5l| and was shown to reduce the scale dependence of the one-loop 



calculation for and B. In Ref. [^, the 2-loop strong and Yukawa corrections 
were shown to have a similar stabilising effect upon the Higgs pole masses. The 
spectrum generator SOFTSUSYl .8.7 used here includes the dominant terms of this 
two-loop subset in both the effective electroweak potential and the calculation of the 
neutral Higgs boson masses. Terms of order hf, hjg^, hlg^, hf, hlhf, h^, h^hl are 
all included. It will clearly be desirable to include other sub-dominant terms from 
Refs. [0, ^ when they become available in a form easy to include in the spectrum 
generator. 

This preliminary investigation sets the stage as concerns the improvements, both 
theoretical and experimental, that are needed if one aims at a prediction of the relic 
density starting from the mSUGRA (high scale) parameters. Otherwise, some of 
these uncertainties should be taken into account as additional systematic errors. 
We then address the question of how precisely, purely within mSUGRA, do we need 
to reconstruct the defining model parameters mo, M1/2, tan/5, to meet the accu- 
racy of WMAP, say. We will give below the detailed procedure of how we arrive 
at the required accuracy on these parameters. This model dependent approach is 
of relevance for analyses at the colliders, especially LHC analyses, where the pa- 



rameters can be extracted from a global fit to a number of observables |5B|, PU[ . In 
some instances, the standard model input parameters such as the top mass or the 
strong coupling are crucial because of RGE effects. We also investigate the required 
accuracy on these parameters. Let us mention in this respect that the uncertainties 



in mft(mb) were investigated in ref. [54| for the funnel region. 

In a second approach, which for short we will call PmSUGRA, we aim at a less model 
dependent approach. The spectrum is still generated within mSUGRA in one of the 
three regions in accord with WMAP. Here however we aim at finding out the accura- 
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cies on the physical parameters at the electroweak scale. This may seem a daunting 
task since the calculation of the relic density involves quite a few channels^ which 
in turn means having access to the entire parameter space of the general MSSM. 
The fact is that, for the three regions under study, the relic density is controlled 
by a quite small set of physical parameters and is rather insensitive to the rest of 
the physical parameters. We foresee this to be the case in most SUSY scenarios 
constrained by WMAP. We therefore start by concentrating on these relevant weak 
scale parameters. To help pin down these parameters we will show, along the slope, 
the relative contribution of those most important cross sections that are necessary to 
predict the relic density. It is educative at this point to characterise how the relative 
importance of the different channels comes about. For the co-annihilation channels, 
the Boltzmann factor is, most of the time, an overwhelming factor. This thermody- 
namical term requires a rather precise knowledge of the mass difference between the 
LSP and the NLSP. Analyses at colliders, in particular the LHC, may have direct 
access to this mass difference rather than on the individual masses. Knowledge of 
the cross sections requires the knowledge of couplings. Some of these are set not 
only by SM couplings but they also involve for example mixing in the neutralino and 
stau sector. In the funnel regions, mass measurements are again crucial due to the 
resonant structure of these contributions. Anyway, in this PmSUGRA approach we 
will, after identifying the relevant parameters, seek to find out the accuracy needed. 

We will generally refrain from addressing the issue of how either the GUT scale 
parameters or the relevant physical parameters in the PmSUGRA approach can be 
extracted in a global fit, or otherwise, from physical observables at the colliders. 
Our main concern is to derive the required accuracy needed to match WMAP (or 
PLANCK). 

3.2 Deriving the accuracies through an iterative procedure 

We ask the question: what fractional change a = \Ap/p\ in an input parameter p 
results in a fractional change r = Auj/uj = 10% for WMAP accuracy? u = Q^oh"^, 
for short. The PLANCK accuracy is obtained by dividing that of WMAP by a factor 
of 5, when the WMAP accuracy is not too large (below 50%, say), so we will usually 
decline to list it for the purpose of brevity. We answer our question by an iterative 
procedure'^. Let us choose a point in parameter space and denote the resulting relic 
density as u^i. Taking an initial guess for oq of 0.1, for instance leads to a certain 
value Uq for the relic density. Taking i = in the equation 

r 

di+l = tti- r- , (3.2) 



^To cover the most general case, micrOMEGAs||55| p3[ includes some 3000 channels or so. 
^Taking the linear approximation and using dlnuj/dhip results in large numerical errors in some 
cases. 
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leads to a prediction for ai. Here, Ui is the relic density that results from using the 
parameter p changed by a fractional amount Oj, so we may iterate Eq, |3]^ up to 
arbitrarily high i. For a given trial in the fractional change of the parameter Oj, the 
fractional factor in Eq. |3.2| is larger than one if Wi comes out to be too small and 
smaller than one otherwise. As i — oo, — a. Thus, we apply Eq. |3.2| successively 
until (oj+i — a.j)/ai < t, where t is the numerical tolerance of the iteration (taken 
here to be 0.001). The resulting value of Oj gives us a good approximation to \Ap/p\. 
We only examine the sensitivity to one parameter at a time and so the rest are 
all set at the values used for the relevant point in parameter space being studied. 
It should therefore be borne in mind that when quoted, precisions on inputs are 
necessary but not always sufficient. For example, the accuracy of several inputs 
must be constrained simultaneously. There is a systematic in the calculation of Qh!^, 
through micrOMEGAs, that we have not taken into account. The error in solving the 
Boltzmann equation could be as high as 1%, not so important when compared to 
WMAP accuracy, but vital when one goes to PLANCK accuracies. This error will 
however be greatly improved by the time we have identified the model as by then we 
will have a tailor made code, rather than a general purpose code, for the calculation 
of the relic density. 

We now turn to the definition of the sensitivity parameters in both the mSUGRA 
and PmSUGRA approach we have just outlined. 



3.3 Analysis in the mSUGRA approach 

We calculate a(mo), a(Mi/2), a(tan/5), a(as(M^)), a(mb(mfc)) by the iterative proce- 
dure described above: given an initial mSUGRA point, we first change e.g. mo, while 
re-calculating the spectrum and couplings. We simultaneously satisfy constraints 
coming from REWSB (Radiative Electroweak Symmetry Breaking) conditions, ex- 
perimental inputs on SM data and the theoretical boundary condition on the high- 
scale SUSY breaking parameters. We denote the new value mfj, which is found by the 
iterative procedure to induce a fractional increase of 10% in u. a{mo) = rriQ/mo — 1 
has analogous expressions for M1/2, tan p. When a is calculated for other parameters, 
mo is set back to its default value. ci(Ao) needs a slightly different definition since 
we will be examining the case Aq = 0. To avoid division by zero in this case, we 
will simply calculate which AAq obtains the 10% fractional increase in u, that is, 
the overall change rather than the fractional one. For reasons detailed later in the 
text, we need to consider a(r^), the total width of the pseudoscalar Higgs boson. 
This is calculated by finding a(mfc(mfc)), which corresponds to a certain a(r^). 
is sensitive to mb(mb) through decays into bottom quarks since the Yukawa coupling 
is proportional to it. So here the use of F^ is traded off for mb{mh)'^. 

^This is done for technical reasons having to do with the way is calculated and used in 
micrOMEGAs. 
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By default, we will consider 10% increases in u; (r = 0.1). There were cases when 
not all of the parameter space considered allowed a 10% increase in uo by varying 
certain parameters. For example, if the required inputs of as{Mz) or mt become too 
large, QCD becomes non-perturbative or the top Yukawa coupling reaches a Landau 
pole before the GUT scale. When this was the case, we chose to examine the value of 
a that results in a 10% decrease in a; (r = —0.1). The sign of r does affect the value 
of a obtained in some cases by a factor of up to 2. Whenever we have used r = —0.1, 
a{p)~ or a{p~) will denote the value of a (ie a superscript "-"). When a(p) > 1 for 
a particular parameter p, we will often choose not to include that parameter in the 
plots (unless it is expected to be very difficult to measure) for brevity. 

3.4 Analysis in the PmSUGRA approach 

Here, we again pick a parameter point in mSUGRA derived from the high-scale SUSY 
breaking terms mo, M1/2, ^0 and which is compatible with REWSB and with the data 
on Standard Model particle masses and gauge couplings. But when we calculate the 
required change a(p) in the parameter p we assume the more general MSSM. As ex- 
plained above, this allows us to leave the strong mSUGRA assumption behind when 
analysing data. The parameter p can be a weak scale parameter that defines the 
effective Lagrangian, examples being the weak scale values of Mi^2, fJ', ^T,nT'ea,eL,f[i,- 
beside the ubiquitous tan /3. These may not be directly extracted from physical ob- 
servables at coUiders, however their advantage is that they set the physical masses 
and the couphngs of the SUSY particles through their mixing matrices. We will, 
however, most often express the accuracies on the physical parameters, such as the 
physical masses. This approach has a direct interpretation in terms of what is mea- 
sured at the colliders. On the other hand in some cases its implementation for 
deriving the accuracies can prove tricky. For example insisting that the neutralino 
mass matrix keeps its tree-level form as is implemented in SOFTSUSY means that one 
can not vary one of the neutralino physical masses while freezing the diagonalising 
matrices that define the couplings. Another example concerns tan (3. For example, in 
the funnel region most of the tan (3 dependence is contained in the Yukawa coupling 
through the pseudo-scalar Higgs exchange. Measuring will get rid of the large 
tan /3 sensitivity. However this will not translate into a measurement of tan /3 since 
F^ depends on other parameters. It is then interesting to ask how sensitive one is to 
the tan f3 which enters the A^^Xi couplings, which might be more difficult to extract 
exp er iment ally. 

4. Co-annihilation 

4.1 Characteristics of the co-annihilation sliver in mSUGRA 

In this region the LSP neutralino is mostly a bino with a mass set essentially by Mi 
up to corrections of order M|/// {/i is large). In the approximation where tan/3 is 
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large, in fact tan /5 > 5 will do, we may write : 



^x?^ ^1-0(^1 + ^) and Ml = 0.417M1/2. (4.1) 
The mass of the right-sleptons in mSUGRA can be approximated by 



2 



rrir 



-ml + 0.152M2/2 ~ + (0.390Mi/2)^ (4.2) 

In the same approximation, the mass of the left-sleptons may be cast into 

m?^ ~ + 0.52M2/2. (4.3) 

Therefore for a small enough mo one can easily understand how m^o ~ ^jj^- More- 
over, with high enough tan/5 and a larger r Yukawa coupling, rrif^ is even lighter 
than the Cr and the fiji. First rrif^ < "mg^^^^ through RGE running effects involving 
third family Yukawa couplings. Second, L — R mixing in the stau mixing matrix 
makes m^^ < m^^. fi is the lightest of the two staus and we define it in terms of 
the current eigenstates with the mixing angle 9f such that fi = cos6'ffR + sin^^f^. 
Approximations for the fi mass and the mixing angle, neglecting the A^- contribu- 
tion, will also prove useful. Using the weak scale parameters, in our scenario we may 
write 



In mSUGRA this can be further simplified as 



tan26'^ ~ ^ ~ — -. (4.5) 

(0.6Mi/2)2 m^o m^o ^ ' 

As a measure of mixing we choose to display our results as a function of cos 2^^, 
since this quantity directly enters the cross sections and decays involving the f's. 
eiJG£;(/^tan/3) is the Yukawa contribution to the RGE running. 

The approximate formulae for the neutralino mass and the fi mass suggest that 
degeneracy occurs for itlq ~ 0.145Mi/2- The mQ— M1/2 constraint imposed by WMAP 
can be made more precise. We can define a slope that parameterises the mo — M1/2 
line where co- annihilation occurs in accord with the WMAP data, see Fig. |I]. We 
will specialise to the case /i > 0, tan /? = 10 and Aq = 0. Following Ref. we take 
a slope "SI" in parameter space: 

mo Ml/2 ^ fMi/2\^ , , 

° - 5.84615 + 0.176374— ^ + 1.97802 X lO^M—^ . (4.6) 



GeV GeV V GeV 
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Figure 1: Slope SI in the mg — plane (full line) and the WMAP allowed region 

(green/light grey). The blue/dark grey area is excluded by direct LEP limit on sparticles, 
the dotted vertical line gives the Higgs mass limit from LEP2 (the region to the left of this 
line is excluded). The red area is ruled out by the cosmological constraint that the LSP is 
neutral. 



This slope has a relic density in rough agreement with the WMAP range. Along 
this slope, Fig. ^ shows the relative contributions of the most important annihilation 
and co-annihilations processes that allow to bring down the relic density to the 
required WMAP level. The physical parameters that enter the most important 
channels will be the most relevant parameters that need to be measured precisely at 
the colliders. The annihilation percentage in the figure corresponds to the proportion 
of neutralinos annihilated by a particular process. We also note that the relative 
importance of some channels changes considerably as the neutralino mass increases. 
It is important to stress that although Xi^i co-annihilation is important, it never 
makes up more than 50% of the overall contribution. For values of AM = m^o — 
iTLf^ ~ lOGeV, it is of the order of the annihilation contribution XiXi ~^ ^ill. The 
largest value of AM along the slope SI displays the onset of the co-annihilation 
region, so one expects the annihilation processes (usually associated to the bulk) 
not to be negligible. As AM gets very small, of the order of a few GeV, they 
are less important than fifi which then dominates. With AM a few GeV, the co- 
annihilation channels involving the smuon and selectron co-annihilations are also 
important. Therefore independent measurement of the smuon and selectron masses 
in model independent analyses will also be required. It is important to note that as 
the mass scale M (set by the LSP mass, say) increases we require smaller and smaller 



- 12 - 




Figure 2: (a) Contribution of the various channels to the relic density in %. (b) AM and 
cos Of as a function of ruf-^ . The abscissa corresponds to M1/2 = 350 — 920 GeV. (The label 
f here stands in fact for fi.) 

AM. This can be understood from the fact that the cross sections which scale as 
oc are less and less efficient and we require new channels to contribute. Past 

M ~ 400GeV this mechanism ceases to be viable. 

Fig. 0(b) shows that AM decreases dramatically from lOGeV to less than IGeV 
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Figure 3: Mass spectrum for some of the SUSY particles along the f co-annihilation slope. 

(that is, less than the r mass) for the higher fi masses. Such f 's will not have a 
two-body decay and might make the reconstruction of some of the needed parameters 
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a quite difficult task. In particular, information extracted from the r polarisation as 



suggested in [£6[ may be nearly impossible to obtain, especially since we are never 



in a very strong r mixing scenario. Indeed, the mixing angle is small throughout the 



slope but at high masses it is tiny. This behaviour can be understood from Eq. O 
for the mixing angle, showing that indeed the mixing decreases with increasing LSP 
mass. We also note that fi ~ '^X3 4 ~ 3m^o on the slope as shown in Fig. ^. This 
figure also helps us to see which processes might be accessible at the colliders. For 
instance at the onset of the co-annihilation region, apart from the production of 
the three lightest sleptons /i/?, e^), there is a window for producing the heavier 
sleptons as well as the associated X1X2 ^ SOOGeV e"*'e~ machine, see also|5^. Of 



course the maximum centre of mass energy will access more channels. It may also be 



more judicious to perform a threshold scan for the hghtest stau as was done in 
however since we also want to access the smuon and selectron parameters a combined 
analysis may prove more appropriate |[40||. For the LHC it is worth mentioning that 



the cascade decay q X2 ~^ '^R ~^ Xi be available which can help in a good 
reconstruction of the model parameters in the context of mSUGRA. 

To close this subsection it is interesting to discuss the coupling Xi^i'^- As seen 
from Fig. this coupling enters in a large number of processes. Defining its right- 
handed part by cr and left-handed part by cl we have, in the hmit of moderate to 
large tan/? we are considering: 

Cr = j= 2 — cosfeifA^ii + smOf — — Nis , 

V2 V Cw Mw J 

3 f ■ a r^w , X . m^tan/5 \ 

Cl = — = sm6'f( — Nu + A^i2) - cosalf — — Nis . (4.7) 

V2 V Cw Mw J 

Nij are the elements of the matrix that diagonalises the neutralino mass matrix, see 
Eq. In the scenario we are in with largish /i, this becomes, at the leading order, 

g sw a , ■ a "^rtan/5 

Cr = -= — 2 cos t^f + sm ijf 



^Cw\ /i 

9 sw ( . n "^rtan/5\ 

Cl = —p^ — smUf — cost^f . (4.8) 

^/2cw \ 1^ J 

Because of the small stau mixing angle and because tan f3 is not too large, the Yukawa 
contribution does not make too much of an impact. We are therefore looking for a 
dependence on the mixing angle coming from tiXi ^ '^7 essentially of the form 
3 cos 26 J. + 5 from the cross sections. 

4.2 Theoretical uncertainties in mSUGRA 

As pointed out earlier, predicting the relic density from the GUT scale parameters 
is subject to theoretical uncertainties such as the order at which the RGE are im- 
plemented. Fig. |a (for M1/2 = 350 - 920 GeV, mo = 73 - 185 GeV) shows how 
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Figure 4: (a) Qh'^ along slope SI. The central line gives the default SOFTSUSY prediction, 
with the broken lines representing the limits due to the scale variation (as defined in 
the text), (b) Effect of approximations along slope SI, the same default MgusY scale 
is assumed here. We show the relic density for (top to bottom): the default calculation 
neglecting SUSY corrections to the tau Yukawa coupling, the default ("full") calculation, 
the default calculation neglecting loop threshold corrections to the neutralino mass, and 
the default calculation neglecting 2-loop terms in the renormalisation group equations of 
gaugino masses. The abscissa corresponds to M1/2 = 350 — 920 GeV. 



the prediction of the relic density changes as one varies the renormalisation scale 
MsusY at each point along the slope SI by a factor 2 in each direction. In passing, 
observe that the slope parameterises the co-annihilation region rather well. Control- 
ling the scale uncertainty will be important in the purely mSUGRA approach since 
these analyses may have to rely heavily on the spectrum evaluator. The band of 
flh'^ values found is shown as the broken lines in Fig. This theory uncertainty 
is roughly ±5% at small , rising to ±15% at higher masses. This effect is due to 
the fact that Msusy varies over a much wider range at higher masses, also that AM 
becomes very small and therefore fractionally more sensitive to Mg^gY■ 

In order to see what in the RGE or the calculation of the pole masses from the 
DR is most important (and thus need to be improved to deal more appropriately 
with the co-annihilation), we look at the effect of some higher order loop corrections 
that affect the neutralino and stau sector. We expect the next order of presently non- 
calculated higher order corrections to have an effect of two orders of magnitude below 
the removed corrections (since they will be suppressed by a loop factor ~ 1/100). We 
concentrate on approximations that have an appreciable effect upon Qh'^. In Fig. 
we show the effect of the important approximations along slope SI by comparing the 
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case when all corrections are included, and when some are omitted. We see from 
the figure that the two-loop terms in the RGEs of the gaugino masses are essential 
to include, as are the one-loop threshold corrections to the neutralino masses. Not 
including these contributions may even make m^^ > m^o and explain why some of 
the curves in Fig. ^ do not extend as far as the curve corresponding to the full 
calculation. Their effect is far greater than the one induced by the SUSY scale. This 
could hint at the necessity of including the three-loop gauge running and two-loop 
threshold effects. Some work in this direction has been performed and the results 
seem to indicate that at least for the LSP mass the three-loop running in the RGE 
has a negligible effect. The threshold loop corrections to the tau Yukawa coupling 
have a small effect. We were unable to study the effect of threshold loop corrections 
to the stau masses, since SOFTSUSY does not yet include them^. According to ref. 
corrections to nif are typically less than 1%. 

4.3 Accuracies of the GUT scale mSUGRA parameters 
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Figure 5: Required fractional accuracies in mSUGRA parameters along slope SI in order 
to achieve WMAP precision. The range of the abscissa corresponds to Mx/2 = 350 — 920 
GeV. The Aq sensitivity is given in absolute terms, AAq, since the reference point is Aq = 0. 

In the mSUGRA scenario, we assume that combined fits from the LHC data agree 
with the mSUGRA predictions, and that mo,Mi/2 and Aq can be constrained from 
observables involving other particles than just those relevant for the co-annihilation 
region. If the theoretical predictions are all under control, this means that we could 
predict the relic density without having accurate information on the stau mass or on 
AM (which might be quite difficult if not impossible). How small should then be 

^We plan to include loop corrections to stau masses in a future version of SOFTSUSY. 
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the precision on the experimentally reconstructed GUT scale parameters to match 
the 10% precision of WMAP? Fig. ^ shows the required accuracy, on both mo and 
Ml/2, to be from around 1 — 3%, with more accurate measurements being required 
for heavier staus. This is not a surprise since M1/2 is directly related to the LSP 
mass and mo contributes directly to m^^ and therefore small changes may make co- 
annihilation not work. More accuracy is needed for higher mass scales since the mass 
difference is very small there and the usual XiXi annihilation is negligible. This is in 
line with the observation that the required precision is essentially set by the Boltzman 
factor, exp —AM/Tf, with Tf ~ m^o/25, as we will discuss later. This also explains 
why the accuracies on mo and M1/2 almost coincide. It rests to see whether such 
an accuracy can be reached considering that from previous analyses mostly done in 
the bulk region, the precision ranged from 1% to 10% depending on the point in 
parameter space, see for example ||53|| . Another LHC analysis [^], also done in the 



bulk region, showed that the precision on the extracted mSUGRA parameters was 
such that a precision of about 3% on the relic density can be achieved. Although 
we believe that, because of the decay chain q —>■ X2 ~* ~^ Xi present also in this 
scenario, some of these analyses can be carried over to the co-annihilation region, 
the small AM value means that some of the accuracies may only be reached near 
the onset of co-annihilation. We also see that Aq should be measured within about 
lOOGeV. The accuracy on Aq should be met in analyses similar to those carried for 
the bulk region, where an accuracy of the same order as we require is found |3^. The 
accuracy on tan (3 is about 9% and is again about what is found for a particular 
bulk region point in a LHC analysis[^]. Sensitivities to the standard model input 
parameters mb, as and rrit are all negligible here. 

4.4 Accuracies on the relevant physical parameters 

The relevant physical parameters here are the masses of all the light sleptons, in par- 
ticular the fi, as well the mass of the neutralino LSP. Once the weak scale parameters 
are set, we take as input parameters for the scalars, mf^,mf2 and the mixing angle 
of the stau cos6'f, as well as the masses of cr and fiR. For the neutrahno one needs 
not only the masses but also their couplings. These involve the diagonalising matrix. 
To be consistent we always start from the set Mi,M2,/i and tan/3, although as we 
discussed earlier the co- annihilation region will be most sensitive to Mi through m^o. 
AM is a crucial parameter and we analyse its accuracy through small changes in the 
physical mass m^^. All other parameters are then fixed. 

Fig. P shows that the mass difference must be measured within just less than 
IGeV. Qualitatively this can be understood from the fact that changing m^^ affects 
those channels where co-annihilation involving fi is important and crucially depend 
on the Boltzman factor. Thus an approximate derivation of the accuracy on the 
relic density can be arrived at by looking at the variation of the Boltzman factor 
dlnQh"^ ~ 6AM /Tf ~ 256AM /m^o ~ 256AM /mr, with small variations caused 
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Figure 6: Required accuracy on the stau-neutralino mass difference obtained by varying 
the fi mass while keeping all other parameters fixed. Accuracy on cos20f is also shown. 
The required accuracy on the mass of Xi is performed by keeping AM constant. 6M = 
6{AM). The abscissa range corresponds to M1/2 = 350 — 920 GeV. 




Figure 7: (a) Keeping AM constant we give the fractional variation in the relic density 
due to a 100% change in fj, and tan (3. (b) Accuracy on the common (right)selectron and 
smuon masses given in terms of S~{AM) = 5{mej^ — ?7i^o) in order to achieve WMAP 
precision. 

by the weight of the other channels. This accuracy can be turned into an accuracy 
on mf^, a(mfj ~ 1 / nif-^^^GeV) , leading to an accuracy that ranges from 0.7% for 
the low ~ ISOGeV to 0.25% for m.^^ ~ 400GeV. cos 26*^ must be measured to 
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within 0.35 (when its value is about ~ 0.92) to about 0.1 (when its value is almost 
1). The reason it is larger at low m^^ is that the stau dependence is almost all con- 
tained in fiXi — > T7, whereas as AM increases, more stau channels are involved, 
for example fifi annihilation through Z-exchange. We now seek the accuracy on 
the overall scale, defined as the accuracy that is required on m^o but leaving AM 
unchanged. This of course means that we also need to change all slepton masses by 
the same change that occurs in m^o. The purpose of the exercise is to determine 
whether all of the sensitivity on the masses comes from AM. Moreover it may well 
be that measuring individual masses, especially the LSP mass, may be more difficult 
in some situations than measuring the mass differences. To determine this accuracy, 
we have chosen to vary the Mi parameter, leaving all other parameters of the neu- 
tralino sector intact. This means however that small changes in both the neutralino 
masses (notably the LSP) and of the couplings of the LSP will be affected. The 
same change in the LSP mass is then applied to all the lightest sleptons such that 
all relevant mass differences and in particular AM, remain unchanged. We find an 
accuracy, corresponding to the 10% WMAP precision, ranging from 20% for the low 
fi mass to 5% for fi masses around 400GeV. This is consistent with a simple scaling 
law for the relic density through the scaling of the cross sections. At the higher end 
of the spectrum, AM is negligible and there is essentially only one mass scale set 
by the neutralino mass. All participating cross sections will scale as 1/M^, once we 
freeze the Boltzman factor. Then the required uncertainty on M = m^o is half the 
uncertainty on the measurement of the relic density. For lower fi masses, terms of 
order AM/M give non negligible effects. Note that this accuracy is about a factor 
20 less precise than the one we obtained by letting AM vary. Therefore the priority 
will be on a very precise measurement of AM although precision measurement of the 
overall scale is also demanded. We can also ask, if AM is measured very precisely 
do we also need to measure tan/5 and n, both of which enter the couplings of the 
sleptons to the LSP? To find any appreciable change in this context requires varying 
these two parameters quite substantially while fixing AM. With a 100% change in 
tan (3, the effect on the relic is below the 2% level, getting smaller at higher masses 
as shown in Fig. |^. Doing the same with /i we obtain a 5% change at low rrifi which 
gets rapidly smaller as we get to higher masses. 

Making sure that we are in the co-annihilation region requires in addition the mea- 
surement of the smuon and selectron masses. Fig. |^ shows that we need to measure 
these to about 1.5% to achieve WMAP precision. This is obtained by varying the 
selectron and smuon mass by the same amount, keeping everything else unchanged. 
Very recent simulations p9|, ^ of the co- annihilation region, done in the context of 
a linear collider, indicate that the relevant parameters AM, m^o,mfj and me^,m^^ 
can all be measured with the required accuracy to meet WMAP precision if one 
is in a region with large enough AM (above 5GeV say). For increasingly smaller 
AM < 5GeV associated with higher fi masses, the situation is more problematic and 
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will need centre-of-mass energies in excess of 600GeV. Future investigations need also 
to study precisely whether the accuracy on the measurement of the mixing angle can 
be achieved in such a scenario. Detailed simulations ||59|| at the LC have been per- 
formed but only in a scenario with quite large AM ~ 80GeV and fi masses around 
ISOGeV, the lower end of slope SI. One can obtain for example cos 26r = 0.987±0.06. 
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Figure 8: Required experimental information along slope SI in the PmSUGRA scenario 
in order to achieve WMAP accuracy on the predicted relic density. Fractional accuracy 
required on fi, Mi. The abscissa range corresponds to M1/2 = 350 — 920 GeV. 



We could also ask what precisions on /i, tan/5 and Mi are needed if we do not 
constrain AM to a fixed value. For Mi, we see in Fig. ^ that we require a few per- 
mil accuracy. This is the accuracy that we found for rrif^ . The latter can be derived 
trivially from AM in Fig. ^ This is another indication that most of the accuracy is set 
by AM. However if one tries to derive the accuracy on /itan/3 without constraining 
the physical masses, by looking at variations on the Lagrangian parameters. Fig. ^ 
shows that one might also need to measure these two parameters quite precisely to 
arrive at a precise determination of the relic density. This shows, when compared to 
our previous analysis, that most of the tan/5 and /i accuracy is contained in AM. 
In fact, the tan/3 accuracy shown in Fig. |^ can be derived from the approximate 
formulae given in Eqs. 4J..4J and 4^. Note in passing that the accuracy on Mi 
is more demanding than that for M1/2 in the pure mSUGRA analysis, although 
we might think at first that since Mi is directly proportional to M1/2, we should 
expect similar accuracies. It should be remembered though that M1/2 contributes 
quite substantially to m^^ through the RGE running in this scenario, see Eq. [4.4| . A 
large cancellation in AM results which explains the weaker dependence on M1/2 in 
mSUGRA. 
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5. Higgs Funnel 



5.1 Characteristics of the funnel region 

Rapid and efficient annihilation can occur through the Higgs resonance. Because 
of the Majorana nature of the neutralino, resonant enhancement is only obtained 
via the pseudoscalar Higgs boson. The heavier CP-even Higgs which may be of 
equivalent mass is completely buried because of the P-wave suppression. In this 
situation A —>■ bb is by far the dominant channel at high tan jS, with some contribution 
from rf. We have checked that for the funnel region we are studying these two 
contributions accounted for more than 98% of the relic density, the rest is annihilation 
into other fermion pairs. This also con&ms that annihilation through Z exchange 
or sfermion exchange is not significant. Therefore we should keep in mind that what 
sets the stage here is the quantity 2M^o — Ma and the width of the pseudoscalar 
since they define the A resonance profile. A precise prediction of the relic density also 
relies on a precise determination of the XiXi^ coupling and also the Abb (and to some 
extent Arf). In fact adding the contributions of the b and r final states amounts, 
to a very good approximation, to expressing these two contributions in terms of the 
total width of the pseudoscalar. To track down the dependence of Ta remember that 
Ta oc M^tan/5^(m^ + m^). The rrif, and tan/3 dependence will be made clearer later. 
As known [|60| the peak structure in the cross section gets smeared once we perform 
thermal averaging in order to arrive at the relic density. Approximations based on 
an expansion in the relative velocity only grossly refiect the result. Keeping this 
important fact in mind, we shall nonetheless, in order to see which most important 
parameters are at stake, set the relative velocity of the neutralino to zero^ . The 
important parameters are then encoded in 

((2m^o)2-Ml)' + riMl 4m.oTA f /M^-2m.o' 
< av >,=o« ^ 4 1 + 1 



(5.1) 



is the total width of the pseudoscalar Higgs while f ^ is its total width into 
non-supersymmetric particles. In our scenario the difference is only due to the con- 
tribution of A — i> XiXij which is negligible as it only occurs very near threshold and 
in any case does not compete with A bb for such high tan/? values. We have 
checked that it constitutes a mere 0.05% of the total width and could thus be safely 
neglected. By default, this SUSY contribution is not added in micrOMEGAs. When 
looking for the accuracy on the total width in PmSUGRA, we will in fact vary the 

^In the situation we will be studying, this approximation should not be too untrustworthy since 
we will be somewhat away from the pole. As we will see (Ma — 2m^o)/rA ranges from 1.4 to 
10. Moreover, the width at tan/? = 50 is not that small. Points directly on the resonance lead to 
annihilation rates that are too low to be within the WMAP range. 
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Abb vertex through a change in rrib so that in effect we are identifying Ta = ^a- In 
the same way, an ag variation in PmSUGRA will affect both Ta and Ta- As we will 
see, the Higgs funnel occurs for rather large values of /i and relatively large Ma- In 
mSUGRA these two parameters are strongly related. The neutralino is still essen- 
tially of a bino nature, but in order to provide a coupling to the A some Higgsino 
mixing is needed. The coupling XiXi^ is, approximately, controlled by 

r(MiS2/3 + /i) , (5.2) 



^?X?X?A OC ^2 _ ^2 ^^^^1^2/3 + /ij 



where the second step of the approximation assumes large tan P and Mi -C /i. 
Ta is an important parameter. Its main contribution is from A bb which is subject 
to large QCD corrections. The latter can be absorbed by the use of an effective b 
quark mass m^-^-^ expressed in terms of the MS b mass at the scale of the Higgs 
boson and some subleading corrections. The Higgs coupling also contains SUSY 
corrections that affect the bottom Yukawa coupling. These can be resummed and 
we parameterise them through the quantity Am;,. In mSUGRA Am^ is also an 
important parameter for the weak scale boundary condition at high tan/5. It is 
derived in SOFTSUSY from 

MMz)WssM = m,iMz)Wi/a + ^ra,), (5.3) 

which re-sums some ~ corrections in the SUSY loop contributions Am^. 
An effective Lagrangian for the couplings of the Higgses can be written as 



eff 

r Ha — 



1 + Arrib 2Mw sin 9w 



^-.,7 cos a / Arribtana 
—Hbb 1 ' 



cos/5 \ tan/5 



^ I . Am-b \ , - sina Am^ 
+z75A66tan/? 1 ^ ] + hbb ' ^ 



tan/52y ' cos/5 \^ tana tan /5y J ' ^ ^ 
where the angle a originates from the diagonalisation of the CP-even Higgs system. 

Anih in the funnel region we are studying is about 25%. Its theoretical deter- 
mination requires the knowledge of a large part of the SUSY spectrum, especially 
the third generation squarks, gluinos and also the Higgsino. All of these may prove 
difficult to measure as these masses are high. One should therefore try to determine 
directly from experiment. As an alternative we may also wonder whether Arrif, 
could be extracted from precision measurements of the light Higgs, h. We feel that 
this will be extremely difficult if one is in the decoupling regime that occurs for large 
Ma > 200GeV. In the decoupling regime, h is essentially SM-like and the Anib 
correction is screened being suppressed by the factor M|/M^ (although large tan/5 
delays this decoupling in hbb, see [^). To wit, normalised to the SM coupling, the 
h and A couplings write respectively, as 



hbb — sin(/5 — a) — cos(a — P)A 



bb 



^ , ^ Arrib \ 7 T tan/5 , ^, 

Abb = Abb[l-^,)^Abb; Abb = ^^^^. (5.5) 
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sin(a — (3) is almost 1 in the decoupling regime and could be determined from the 
measurement of the ZZh coupling, for example. It would be interesting to examine in 
such a scenario whether one could combine the extraction of Ahi, through a precision 
measurement oi h ^ bb at the linear collider and the LHC measurements of Ma 
in order to determine A ^ bb more precisely. There is also an equivalent of Arrib 
for A ^ Tf, however it is much smaller since the corresponding Yukawa coupling is 
much smaller and because it does not feel the important QCD correction. 

In our discussion of the Higgs funnel we will take tan (3 = 50, Ao = 0, fi > 0. We 
parameterise the funnel region through the slope S2, defined as a cubic, see Fig. ^a^: 




where q = {814.88, -2.20022, 3.30904 x 10'^, -1.05066 x 10"^}. 




Figure 9: a) Slope S2 in the ttt-q- M]^/2 plane (dashed line) and the WMAP allowed region 
(green/light grey). The red area is ruled out by the cosmological constraint that the LSP 
is neutral, b) Mass spectrum of the relevant particles in the funnel region. M1/2 is in the 
range 250-llOOGeV. 

An important remark is that the funnel region has a rather heavy spectrum, see 
Fig. ^ In particular the sfermions are very heavy making it impossible to produce 
any sleptons at a 500GeV machine. More energy is needed for this. There is some 
chance for the associated production of X2 (through X1X2 production), which could 
serve as a good measurement of m^o and indirectly of the selectron mass. On the 
other hand, the 77 option of a linear collider could bring important constraints on the 

^We found that a quadratic was not sufficient to contain flh^ within the 2cr range. 
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mass and couplings of the A as well as its width. Combined with a determination 
of II from the LHC, one could reconstruct the parameter space that defines the 
funnel region through m^o (observe that m^o tracks fi very well). The LHC can 
also measure the mass of A, possibly up to ITeV, and might probe Ta through 
A /i/i. For future discussion it is interesting to realise that Ta/Ma here is about 
3%, slightly larger than the corresponding quantity for the Z lineshape. Note also 
that {Ma - 2m^o)/TA ~ 1.4 - 10. 

5.2 Theoretical uncertainties in mSUGRA 

The funnel region involves a heavy spectrum. Since it may not be possible to have ac- 
cess to all the relevant parameters, one may have to partially rely on some mSUGRA 
assumptions. It is therefore important to inquire about the theoretical uncertainties 
in mSUGRA when predicting the relic density. 

Fig. p!0| a shows the scale dependence of Qh"^ along the funnel slope S2 which, in 
passing, approximates the WMAP data rather well. The scale dependence increases 
from 5% at low values of M^o to 20% at higher values. Mg^gY ~ {1-5 — '2}Msusy 
did not yield physical points, possessing unstable saddle-point electroweak vacua 
(M^^M'gjjgy) < ^) ' Such values of Mg^gY were therefore not included in the scale 
variation. Fig. p!0| b shows the dramatic effects, especially for increasing Ma, of some 
higher order loop contributions in the spectrum calculation for slope S2 and a fixed 
value of mfc(mfe). Not re-summing the Arub terms has a drastic effect. Taking into 
account only the one-loop terms (instead of two-loop terms in the default calculation) 
for the scalar Higgs ^ , m| or gaugino masses in the RGE is not sufficient at all. 
For ITeV pseudo-scalar mass, the 2-loop QCD correction to the top mass is essential. 
Some of the effects of the approximations can be counter-balanced by changing the 
mb{mb) input. Fig. p!0|c. 

To see if most of the dependence can be reduced if some physical observables are 
fitted, we take the variable {Ma — 2M^o)/rA that describes the A resonance profile. 
The illustration is made for a single point on S2: mo = 1000 GeV, M1/2 = 1100 GeV 
(corresponding to the extreme right hand side of Fig. p!0|a). Fig. p!0| d shows that, 
as one varies the scale, there is a peak in the relic density at the same location as 
in the variable {Ma — 2m^o)/r^ describing the A resonance. This indicates that a 
measurement of this ratio could help pin down the uncertainties. 

5.3 Required accuracies on the mSUGRA parameters 

In Fig. [TT|a, we see the accuracies required on measurements assuming mSUGRA is 
correct, rrit must be measured at the level of 0.6-0.2 GeV, depending upon position 
along the slope. as{Mz), mb and tan/5 must be known at the per-mille level. Most 
of this dependence is from the extraction of the b Yukawa coupling. As we have 
seen this could be controlled by measuring Ma, ^a- To match WMAP accuracy a 
2-5% measurement of rriQ is required, while M1/2 needs to be known at the 1% level 
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Figure 10: Relic density along slope S2. (a) The central line gives the default SOFTSUSY 
prediction, with the broken lines representing the limits due to scale uncertainty. The 
range of the abscissa corresponds to M1/2 = 250 — 1100 GeV. miy{mh) = 4.23GeV. (b) 
Effect of different approximations in the RGE, see text for details, (c) The dependence 
on the input value mh{mh) including treatment of higher orders related to the b mass 
(and coupling Yukawa) for = llOOGeV on mSUGRA slope S2. The meaning of the 

different curves is explained in the text.(d) Scale dependence of the relic density compared 
to the scale dependence of the relevant parameter controlling the dominant contribution 



of the A resonance for = llOOGeV GeV on slope S2. x is defined in Eq. 3.1. 



(although there is a region where an 0{1) measurement of M1/2 would suffice). Aq 
must be measured to 100 GeV accuracy. We have also attempted to quantify how 
the accuracy needed on mb(mb) from WMAP translates into a precision on Mj^. 
For this exercise we of course still need to assume that all other SM and mSUGRA 
input parameters are known precisely so that a theoretical prediction of Ma within 
mSUGRA is possible. We find that one needs to measure Ma to less than 5GeV. 
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Figure 11: (a) Accuracies a required to achieve WMAP precision along slope S2 in the 
mSUGRA scenario, (b) Pseudoscalar Higgs width Ta and required precision upon in 
STa along slope S2 in the mSUGRA scenario. The range of the abscissa corresponds to 
Ml/2 = 250 — 1100 GeV for both (a) and (b). This plot applies to the accuracy needed 
on mi) when all other SM and mSUGRA input parameters are set (or known precisely). 
However we translate the mh{mh) accuracy into the precision needed on Ma- 



This would presumably require the measurement of a line-shape by (for example) 
A IQ. If that were possible we could then use Ma as a trade-off for 

mb{mb) in mSUGRA. On the other hand we found that the width is insensitive 
to mb{mb) in mSUGRA. This should not be taken to mean that with measurements 
of Ta and Ma alone we can predict the relic density within WMAP precision. 



5.4 Pmsugra 

As we have already observed, the relic density in the funnel region derives essen- 
tially from XiXi ~^ {bb + tt) through pseudo-scalar Higgs exchange. The important 
physical parameters are Ma, T^, m^o as well as the parameters that enter the XiXi^ 
coupling which are mainly controlled by fi, see Eq. p.l| . For all purposes yU can be 
equated with m^o, see Fig. |^. Let us therefore start by finding out the accuracy 
needed on this physical parameter first while keeping all other parameters fixed. To 
arrive at this accuracy we simply change the value of /x at the Lagrangian level given 
at the weak scale. Strictly speaking this also changes the value of the LSP mass. In 
the situation we are in, Mi,Mz -C /i,tan/3 ^ 1, the induced change in m^o will be 
a very small change as can be seen from Eq. ^]l]. More important, a change of will 



impact directly on the LSP coupling to the pseudoscalar Higgs. Fig. corroborates 
this observation. We see that one needs an accuracy of about 5% on fi for a 10% 
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Figure 12: Accuracies a required to aciiieve WMAP precision along slope S2 in the 
PmSUGRA scenario. The range of the abscissa corresponds to = 250 — llOOGeV. 

2m^o — Ma is also shown on the right hand axis. 

A.1 



WMAP precision, which corresponds to the fact that Qh"^ oc /x^. The evaluation of 
the accuracy on F^, independently of any of the other relevant parameters, is trig- 
gered through a change in mj,. We find an accuracy on the total width which is of 
the same order as the WMAP accuracy, in accordance with Eq. |57^. We have also 
looked at the accuracy needed on as{Mz), assuming one knows the other parameters 
entering the calculation of Ta- We confirm that the accuracy needed is proportional, 
and of the same order, to that found for P^. The accuracy on the LSP mass is arrived 
at by varying the bino mass Mi, although this also very slightly changes the XiXi^ 
coupling. We find an accuracy on m^o ranging from 2% for m^o ~ lOOGeV to 0.2% 
for m^o ~ 450GeV. This can also be translated into an accuracy on the A profile 
parameter (2m^o — M4) which we find to be, as expected, of the same order as fi: 
5% for WMAP accuracy. The accuracy of Ma is derived by varying Ma but note 
that this also changes Pa, oc Ma- As expected from Eq. ^]T] we get an accuracy 
on Ma which is practically the same as the one on m^o, from 2% to 0.2% for WMAP 
precision. Experimentally, this will be quite demanding, especially for the largest 
of pseudo-scalar masses allowed in our funnel region. It is important to add that 
once we choose the set M^, P^, m^o, /x(m^o), the tan/3 dependence is very mild. To 
confirm this we have tried to "decouple" the tan j3 dependence in the vertex Aff 
from the one in the coupling of the neutralino. We do this by taking a point with the 
default value of tan j3. We then change tan (3 and record the changes in the neutralino 
mass matrix. This new neutralino mass matrix is then substituted in the original 
default point, leaving all other variables unchanged. We find for example that a 10% 
change in tan/5 does not affect the relic density beyond 1%. Although negligible for 
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WMAP accuracy, requiring PLANCK accuracy means that a good determination of 
tan P might still be needed. 

6. Focus point 

6.1 The focus point landscape 




Figure 13: a) WMAP allowed region in the tjiq — M1/2 plane (green/light grey) and slope 
S3 (dashed line), b) Mass spectrum in the SOFTSUSY focus point region along slope S3, 
tan P = 50. 



The focus point region ||6^ , ^ corresponds to high values of mo near the bound- 
ary of viable electroweak symmetry breaking (see Fig. ^a) where the value of n 
drops rapidly. When fi ~ Mi,M2, the LSP has a significant Higgsino fraction, fur- 
thermore the next-to-lightest sparticles (^2 or x^) also have a significant Higgsino 
component and are not much heavier than the LSP as shown in Fig. p!B|b. Thus 
co-annihilation channels are favoured. However, co-annihilation should not be too 
efficient, otherwise the relic density is less than what is measured. The currently 
acceptable focus point region obtained with SOFTSUSY, though requiring fi smaller 
than in the co- annihilation and funnel regions, does not lead n < Mi, M2 but rather 
Ml < fi < M2. This still means that the LSP has a sufficient Higgsino component, 
of about 25%, so that stronger couplings to the Z and W take place, in particular 
through their Goldstone component, and to some extent through Higgs bosons. 

We will examine a high tan/3 = 50 focus point compatible with WMAP as 
obtained with SOFTSUSY and take a slope "S3" in parameter space: tan P = 50, ^ > 



- 28 - 



0, Ao = and 

mo Ml/2 ,/Mi/2\^ , , 

—^7 = 3019.85 + 2.6928-^ - 1.01648 X lO^M-^ . 6.1 
GeV GeV \ GeV / 

with Ml/2 in the range 440 — lOOOGeV, see Fig. ITla. This means that all sfermions 
are far too heavy (in excess of ATeV or so) to be accessible at any of the planned 
colliders. The interesting feature of the spectrum is shown in Fig. |l3|b. The LSP mass 
ranges from about 150 to 350GeV with the Higgsino-like neutralinos being some 100 
to 50GeV heavier. Their mass differences decrease for the higher masses and thus 
the Higgsino content of the LSP increases. The smaller mass difference suggests that 
co-annihilation may be more important for higher scales. The wino-like gauginos are 
above 400GeV. From the point of view of a linear collider, an energy in excess of 
800GeV is needed to unravel some properties of this scenario. Although not shown 
on the plot, let us mention that the pseudo-scalar has a mass in excess of ITeV and 
will most probably not be directly probed at the LHC. The gluino mass is also above 
ITeV but, unlike the pseudo-scalar Higgs, does not enter directly in the prediction 
of the relic density. 
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Figure 14: Contribution of the various channels in the focus region. "Annihilation" 
stands for all the channels which we give individually as tt,bb,W~^W~ , ZZ, Zh and hh. 
"Co-annihilation" stands for X1X2' X1X3, Xixf • 

The contribution of the most important channels to the relic density along slope 



S3 are displayed in Fig. IHl First of all, note that in our scenario annihilation into 
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tops is always the largest contribution. This comes about through the Goldstone 
component of Z exchange, which picks up the large top Yukawa coupling. Because 
of the relatively large Higgsino component, the neutral Goldstone boson couples like 
MzMxj (Ml — /i^), with /z not much larger than M\. We also note that we get 
a non negligible contribution from annihilation into 66. This contribution is due 
to A exchange which, contrary to the tt contribution via Z exchange, gets a large 
tan /? enhancement. Properties of this contribution are similar to those in the funnel 
region, though here it is diluted by comparison because of the much larger pseudo- 
scalar masses involved and also because we are far away from the pole. Although 
the XiXi^ coupling is slightly larger than the XiX^Z coupling, there is a factor of 
yu/Mi between the two couplings, annihilation through A exchange is suppressed 
by the large mass of the Higgs in the propagator, compared to the mass of the 
Z for the ti channel. Note the other annihilation channels into WW, ZZ, Zh, hh 
which account for 10%. For LSP masses around 350GeV, co-annihilation with the 
heavier neutralino starts to contribute since the Mi,/i mass difference gets smaller. 
From these observations, the fact that quite a few channels and mechanisms are 
contributing and that the masses involved are quite high will make a precise probe 
of this scenario rather difficult for the colliders. 



6.2 Theoretical uncertainties 

Theoretically, the focus point regime is also fraught with difficulties. The position 
of the focus point region is extremely sensitive to the value of the DR MSSM top 
Yukawa coupling ht |6^, 0, which differs for different RGE codes |5^. ht is fixed by 



the MSSM value of the running top quark mass. In order to obtain this, supersym- 



metric and Standard Model corrections are subtracted from the top pole mass [4£ 
The way in which this is done differs in terms which are formally of higher order 
in the various codes. For example, some codes use pole masses for particles in the 
loops whereas some use running masses. The radiative corrections are sometimes 
calculated at different scales. For the case of SOFTSUSY, values of the top quark mass 
lower than the present la range are necessary to reach the region where the relic 
density is in agreement with the WMAP constraint if tan (3 is intermediate (e.g. 10- 
30). This is the reason we concentrated on a high tan/? = 50. We have checked that 
the numerical results are similar (along a different slope consistent with the WMAP 
measurement) for tan/? = 10 and rrit = 172 GeV. 

It can be seen from Fig. |l5|(a) that the prediction of the relic density for the de- 
fault SOFTSUSY prediction (the unbroken line) is consistent with the WMAP measure- 
ment along S3 between M1/2 = 440 — 1000 GeV (corresponding to mo = 4185 — 5611 
GeV). However, when Msusy is varied by a factor of 2 in each direction, the predic- 
tion can increase by up to 20% and decrease by up to 30%. This illustrates that the 
theoretical uncertainties are large for the focus point regime. 
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Figure 15: {a)Qh'^ along slope S3. The central line gives the default SOFTSUSY prediction, 
with the broken lines representing the limits due to scale uncertainty. The range of the 
abscissa corresponds to M1/2 = 440 — lOOOGeV. (b)Relic density dependence upon different 
approximations in SOFTSUSY for M1/2 = 495GeV, Aq = 0, tan/3 = 50 and > 0. The 
curves refer to the following approximations: (full calculation): the default full calcula- 
tion, (2 loop QCD): the default calculation with the 2-loop QCD corrections removed from 
the extraction of "itA/^sA/ (n-t), (EW/Higgs): the default calculation with the one-loop 
electroweak and Higgs corrections removed from the extraction of nT'tf/ssAii'^t)^ (neutrali- 
nos): the default calculation with the one- loop neutralino corrections removed from the 
extraction oimt^j^gj^,^{mt), (charginos): the default calculation with the one-loop chargino 
corrections removed from the extraction of mt^^^jyj(mf ), (2-loops Higgs RGE): the default 
calculation with only 1-loop RGEs for the Higgs potential soft SUSY breaking parameters 
rn^j^, m^^, m|, the 2-loops Higgs RGE are thus removed, (non-resummed mh)'- the default 
calculation with no re-summation in SUSY mf, corrections. 

We now examine the theoretical uncertainty on the predictions in more detail. 
We use various approximations when calculating the sparticle mass spectrum and 
plot different values predicted for VLh"^ as a function of mo, for M1/2 = 495GeV, 
Aq = 0, tan/3 = 50, /i > 0. We are therefore moving away from slope S3 to show 
that different approximations do lead to different parameterisations. The results are 
displayed in Fig. |l^(b). We can see the extreme sensitivity to loop effects from the 
vastly different values of tuq that agree with the WMAP constraints from the figure. 
The Higgs mass parameter RG evolution becomes very steep near Msusy, and so the 
sensitivity to the number of loops used to evolve is high. The high sensitivity to the 
extracted value of ht is shown by the high sensitivity to the approximations in the 
calculation of ^tA,issMi''^t)- At high^ tan/5 = 50, hh, the bottom Yukawa coupling 
^For the tan/3 = 10 slope, there is no visible difference between the calculation using re-summed 
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is large and affects tlie running of m|^^ appreciably. Its extraction, as we liave seen 
in tlie previous section, requires knowledge of Am;,. Tliere was no visible difference 
from the default calculation by using only one loop running for the gaugino masses 
Ml 2,3 instead of the default two-loop running, or removing all 2-loop terms while 
calculating the Higgs masses and electroweak symmetry breaking conditions. 

6.3 Accuracies in mSUGRA 
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Figure 16: Fractional accuracy required to achieve WMAP precision in the mSUGRA 
scenario along slope S3: (a) Standard Model inputs, (b) mSUGRA inputs. The range of 
the abscissa corresponds to M1/2 = 440 — 1000 GeV. 



We now turn to the issue of sensitivity of the prediction to the inputs in the 
focus point region. As explained in section 2, we imagine that fits to LHC data have 
identified the mSUGRA SUSY breaking parameters. This would then allow one to 
derive the entire set of the weak scale parameters, in particular all of those which 
are of relevance for the calculation of the relic density. As usual it is not only the 
mSUGRA SUSY parameters at high scales which are needed but also the SM input 
parameters. Considering what we have just seen concerning the crucial role played by 
the top Yukawa coupling both in the RGE and in the ti annihilation, we expect that a 
very precise input top pole mass is required as weh as as{Mz), one of the parameters 
that strongly affects the extraction of ht^^^j^^ from nit. Fig. |^a shows that there is 
indeed a huge sensitivity to rrit and even ^^(M^). a{mh{mj,)) ~ (9(1 — 15)% is shown 
as an insert in the figure^. Values of a{'mt) ~ 10~^ mean that empirical uncertainties 

or non-re-summed SUSY m?, corrections. 

^For a tan/3 = 10, mt — 172 GeV focus point sfope, a{mb{mb)) > 1.0. 
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of around 20 MeV on rrit would be required for WMAP precision. Studies suggest 
that S{mt) ~ 30 MeV is possible for the statistical error at a suitable linear 



collider, but that theoretical errors of order lOOMeV are likely. as{Mz)s^M must be 
measured to better than an order of magnitude more precisely than current data 
allow. 

Fig. shows that rriQ must be measured to a fractional precision of better than 
0.2% for WMAP accuracy. Such an accuracy is inconceivable at the LHC because, 
among other reasons, the production cross-sections for the multi-TeV scalars present 
in the focus-point region are tiny. In mSUGRA, mo determines (among others things) 
the value of which determines the amount of Higgsino content. M1/2 also affects 
the determination of as well as the mass of the LSP. We find that M1/2 should 
be measured to better than 0.5% for WMAP accuracy, and Aq to better than 180 
GeV, provided rrit is not quasi-degenerate with M^o. The sudden deterioration in 
the accuracy for M^o ~ rrit is due to the approach to the t threshold. A similar 
behaviour is seen in the required accuracy of tan/?. The measurement of tan/3 in 
mSUGRA will be very demanding, due not only to the importance of the tt vertex 
but also the extraction of the top Yukawa that enters the RGE. 

6.4 Accuracies in PmSUGRA 
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Figure 17: Required fractional accuracy upon MSSM parameters in the PmSUGRA sce- 
nario along slope S3. The abscissas correspond to M1/2 = 440 — 1000 GeV. 
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Fig. 17 displays the required accuracies a for several parameters in the Pm- 
SUGRA scenario. The relevant parameters are the weak scale values of the elements 
of the neutralino mass matrix. In particular, Mi and are important since, as dis- 
cussed already, they determine the Higgsino component of the LSP and how large 
the couplings of the LSP to the Higgs bosons, W and Z are. tan/3 enters in the 
contribution of the hb annihilation and in the neutralino couplings. Since M2 is not 
completely decoupled its effect is expected to feed in as well. All of these parameters 
are essential for the co- annihilation processes that set in for higher LSP masses. Fi- 
nally, there is also dependence upon Ma- Before we discuss these accuracies and in 
order to highlight the contrast with the mSUGRA approach, we comment on the ac- 
curacy needed on mt. mt enters through the coupling of the neutral Goldstone boson 
as (taking the dominant contribution) when the top threshold sets in. We find 
a sensitivity to rrit which is typically of the same order as the precision on the relic 
density. The present empirical precision on is sufficient even for the PLANCK 
benchmark. At the top threshold one requires a much better accuracy. On the other 
hand, the accuracy is less demanding for higher LSP masses since the contribution 
of the ti annihilation is smaller. A rough estimate on a{mt), away from threshold, 
is a{mt) ~ 1/(2q) in units of the precision on the relic and where q is the relative 
contribution of the ti channel, q can be read off from Fig. ^ Similarly in the whole 
range of M1/2 considered, a{mb) > 3.0 meaning that even the present empitical rrib 
errors reach PLANCK accuracy. The sensitivities of /i. Mi require fractional preci- 
sions of 1% for WMAP accuracy. This rather demanding accuracy originates more 
from the couplings of the LSP to the Goldstone and A pseudo-scalar bosons. As we 
show in the figure, the accuracies on Mi and fi can be converted into accuracies on 
m^o and ttl-^o. The accuracy needed on M2 is an order of magnitude worse, though it 
is still non trivial to obtain empirically. Therefore we see that if one can reconstruct 
the neutralino mass matrix, which we think can be done rather precisely at the LC 
provided there is enough energy to produce all of the neutralinos, this scenario can 
be very much constrained. It could be interesting to study whether this could be 
done without direct production of the heaviest (wino-like) neutralino and how a com- 
bination of LHC and LC data could be performed. A slight problem still remains as 
concerns the M^ dependence and its required accuracy. In this scenario Ma is very 
large, in excess of ITeV and could be difficult to access at the LHC or a sub-TeV 
LC. We find that at the lower end of the focus region, a change of 50% or higher in 
Ma only amounts to a maximum 5% change in the relic density. However, at the 
higher end, though Ma ~ l.5TeV, a 30% change in Ma affects the relic density by 
10%. The reason that a larger (positive) change does not deteriorate the precision 
is due to the propagator suppression of Ma which means that the A exchange no 
longer contributes for too high values of Ma- 
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7. Summary and Conclusions 



Assuming mSUGRA, we have examined the requirements for measurements of the 
MSSM spectrum in order to make rehable predictions of the rehc density Qh'^. If 
such requirements could be met, the assumptions that go into the derivation of Qh^ 
may be tested. 

We find that theoretical uncertainties in the mSUGRA prediction of Qh'^ assum- 
ing standard cosmology are too high compared with the WMAP prediction. In the 
co-annihilation scenario, the theoretical uncertainty in Qh!^ (estimated by scale de- 
pendence) is of the same size to twice as big as the WMAP accuracy. The situation 
is worse in the problematic focus point and higher mass Higgs funnel regions, where 
the theory scale variation is around five times bigger than the WMAP accuracy. In 
the lighter Higgs funnel region, the scale variation is smaller than the WMAP ac- 
curacy. The theoretical uncertainties are not a serious problem, since higher order 
calculations should be possible given enough time, effort and motivation. If SUSY 
particles are discovered at a future collider, it seems highly likely that this task will 
be accomplished. 

The experimental information necessary in order to achieve a precise prediction 
of the relic density is demanding. If one assumes mSUGRA, fits must constrain 
Ml/2 and mo at the one or two percent level and Aq to roughly 100 GeV, for the 



co-annihilation region for example. It is estimated in a bulk-region test case [14 



that at the LHC this accuracy is not quite achievable, even when the parameters are 
quite favourable (for example mo = 103 ± 8 GeV, M1/2 = 240.0 ± 3 GeV at post- 
WMAP benchmark point B), but is an order of magnitude too large at other points 
(mo = 400 ± 100 GeV and M1/2 = 400 ± 8 GeV at LHC points 1,2 fUj). Combining 
complementary information from the LHC and a future linear collider should provide 
the necessary precision, however |]64| , [65| , |66| . Not only must one experimentally con- 
trol the soft SUSY breaking parameters assuming mSUGRA, but Standard Model 
inputs must also be controlled if the calculation of the MSSM spectrum is to be 
trusted. In the co-annihilation region this is not a big problem, since the parts of the 
spectrum relevant to the prediction of Qh'^ (the lightest neutralino and stau) do not 
have a large sensitivity to the Standard Model inputs. The sensitivity in the focus 
point region to the top Yukawa coupling hf is so high that m^ must be measured to 
better than 20 MeV. This is smaller than the strong scale of the QCD interaction 
and therefore cannot be perturbatively defined with such precision, rendering the 
approach practically hopeless. In the Higgs funnel region, m^ must be measured 
to around 200 MeV, which may be achieved at a future linear collider facility, but 
mb{mb) must be constrained to the 2-10 per mille level because of its effect on the 
mass of the pseudo-scalar Higgs, which participates in the s-channel to neutralino 
annihilation. This will not be possible since it is much smaller than non-perturbative 
effects, but by measuring the mass of the A at the percent level, the sensitive de- 
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pendence upon mb{mb) can be removed. For the heavier decoupled values of m^, a 
linear collider with sufficient centre of mass energy will be required since the LHC 
will not provide enough precision on measurements of the A bosons [p7[ . 

Rather than relying upon mSUGRA to make predictions for the spectrum, a 
more general strategy is to first determine what regime one is in. If quasi-degenerate 
lightest stau and neutralinos are discovered for example, it would indicate a co- 
annihilation regime. In that case the mass splitting (less than 12 GeV) must be 
measured to better than IGeV in order to achieve WMAP precision upon Qh"^ for 
LSP and fi weighing between 150 to 450GeV. In the focus point regime (indicated 
by the measurement of gauginos but non-observation of scalars at the LHC), the 
extreme sensitivity to rrit vastly reduces if one no longer has to assume mSUGRA 
and relies on individual mass measurements instead. However, a measurement of /i 
and Ml to 1% accuracy is required and will require copious production of neutralinos 
and charginos at a linear collider with sufficient centre of mass energy |68|, |69| since 



the masses involved are high. In the Higgs funnel scenario, the quantity 2M^o — Ma 
must be measured very precisely requiring a measurement of Ma and m^o with a 
precision ranging from 2% for LSP masses around lOOGeV to 0.2% for LSP masses 
around 450GeV. 

The analysis we have performed in this paper was done in the mSUGRA frame- 
work. WMAP now allows only three scenarios in mSUGRA which are quite contrived 
and are challenging both from the theoretical and experimental point of view. Going 
beyond mSUGRA, for example in a scenario with non universal gaugino masses and 
in particular a wino-like LSP, would be subject to less theory uncertainty and might 
call for accuracies that are less demanding for the colliders. It is also important 
to stress that we have assumed that the standard derivation of the relic density is 
correct. However, one should not forget that most alternatives that affect the cosmo- 
logical assumptions usually predict orders of magnitude differences with the orthodox 
approach, far beyond the few percent precision we have taken as a benchmark. This 
means that perhaps some of these alternatives could be dismissed without requir- 
ing too much precision. On the other hand it would be interesting, assuming such 
schemes to be correct, to inquire about the level of accuracy on the particle physics 
parameters that the colliders must meet to corroborate some of these non-standard 
cosmologies. 
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A. Inputs 



Except where explicitly mentioned, we use the following default Standard Model 
parameters as inputs |j70|: 



as{Mz)f^i = 0.1172, aiMz)f^j = 1/(127.934), m,im,f'^ = 4.23 GeV, 

TTir = 1.7777 GeV, Mz = 91.1876 GeV, rrit = 175 GeV. (A.l) 

A subscript SM denotes the fact that the quantity has been derived assuming Stan- 
dard Model field content. mb{mb)^^ is derived from experiment in 5-flavour QCD. A 
MS superscript denotes a quantity calculated in the modified minimal subscription 
renormalisation scheme. The weak mixing angle is derived from the muon decay 
constant = 1.16637 x 10^^ GeV'^. 

The mSUGRA parameters are specified as tan/?, the ratio of the two Higgs 
vacuum expectation values, the universal gaugino mass M1/2, the universal trilinear 
scalar coupling Aq and the universal scalar mass mg. Universality is imposed upon 
the soft SUSY breaking terms at a scale Mqut which is defined by 

giiMcuT) = 92iMGUT), (A.2) 



gi and g2 being the electroweak gauge couplings in the modified DR scheme, gi is 
related to the Standard Model hypercharge coupling g' by the GUT normalisation 
gi = a/5/3(7'. g3 is set by the as{Mz)s^M input and is not unified with gi 2- We will 
take the Higgs potential parameter /i > 0, with its magnitude being fixed by the 



electroweak symmetry breaking conditions [32 
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